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Thomas  F.  George 
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Abstract 

Two  features  of  light  scattering  from  a  deep  lossless  metallic 

sinusoidal  grating  are  considered  in  the  limit  g/d  -*  •  ,  where  g  and  d 

are  the  height  and  periodicity  of  the  grating,  respectively.  It  is 

found  that  the  surface  plasmon  dispersion  relation  is  comprised  of  two 

flat  bands  with  a  frequency  gap  4 u/u>  -  C/2  -  1)/2,  where  w  is  the 

P  P 

volume  plasmon  frequency.  For  the  local  field  enhancement 
distribution  at  the  grating's  surface,  the  reults  show  the  existence 
of  two  qualitatively  distinct  domains,  i.e.,  A  >>  d  and  A  <<  d,  where 
A  denotes  the  radiation  wavelength.  In  both  domains,  however,  the 
local  field  enhancement  i3  larger  at  the  bottom  of  the  troughs  than  at 
the  top  of  the  peaks.  The  dressed  Rayleigh  expansion  is  used 
throughout  for  the  analysis. 
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1 .  Introduction 

Light  scattering  from  a  deep  grating,  i.e.,  a  grating  for  which 

6  -  2irg/d  >>  1  where  g  and  d  are  the  maximum  (and  minimum)  extension 

and  periodicity  of  the  grating,  respectively  (see  Fig.  1),  is  a 

relatively  newly  studied  physical  phenomenon.  Moderately  deep 

gratings  can  now  be  fabricated,  e.g.,  by  photodissociation  of 

1  2 

organometallic  molecules  at  a  substrate.  ’  The  extremely  deep 

grating  case  (6  ■*  «) ,  which  is  at  the  focus  of  this  work,  should  be 

considered  as  an  idealized  limit.  This  limit  is  complementary  to  the 

shallow  grating  limit  (6  <<  1),  which  is  realized  in  many  physical 

systems  and  has  been  extensively  explored  over  the  last  century. 

The  extremely  deep  grating  regime  is  qualitatively  different 

3 

from  the  shallow  grating  regime.  It  can  be  characterized  as  a  strong 

coupling  situation  in  the  sense  that  the  many  possible  Bragg 

reflections  (Fig.  1)  interfere  strongly  with  each  other,  giving  rise 

to  a  scattered  field  of  new  qualitative  features.  By  contrast,  the 

shallow  grating  regime  can  be  analyzed  in  terms  of  a  few  'weak'  (in 

comparison  with  the  specular  reflected  light)  Bragg  reflections. 

Unfortunately,  attempts  to  solve  the  Maxwell  equations  for  deep 

gratings  using  the  so-called  Rayleigh  expansion,  which  works 

beautifully  when  6  <<  1  ,  have  encountered  numerical  instabilities. 

These  difficulties  have  motivated  the  introduction  of  several 

1  5  6 

alternative  schemes,  ’  ’  all  of  which,  however,  share  difficulties  in 

3 

the  6  >>>  1  regime.  In  another  paper  we  have  suggested  yet  another 
scheme  --  dubbed  as  the  dressed  Rayleigh  expansion  —  which  overcomes 

3 

these  difficulties.  We  have  also  explicitly  demonstrated  that  this 


j 
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scheme  converges  for  a  sinusoidal  grating  (SG)  of  arbitrary  Here 
we  apply  the  dressed  Rayleigh  expansion  to  study  two  aspects  of  light 
scattering  from  a  very  deep  metallic  SG  grating,  namely,  the  surface 
plasmon  (SP)  dispersion  relation  and  the  distribution  of  the  local 
field  enhancement  along  the  grating's  surface. 

Our  results  are  valid  in  the  3  -»  ®  limit  and  under  several 

approximations  which  are  spelled  out  as  we  go  along.  For  the  SP 

dispersion  relation  we  find  two  flat  bands  (see  Fig.  2).  The 

/2-1 

frequency  band  gap  Au  is  given  by  Aui/ui^  *  — —  »  .207,  where  is  the 
volume-plasma  frequency  of  the  metal.  We  interpret  the  two  bands  to 
reflect  two  classes  of  localized  surface  modes  —  those  for  which  the 
field  maxima  are  in  the  peaks  and  in  the  troughs  of  the  grating 
respectively.  A  similar  situation  occurs  in  standard  band  structure 
studies.  The  fact  that  we  obtain  dispersion  bands  rather  than  the  usual 


dispersion  curves  for  3  <<  1  ia  in  keeping  with  studies  of  moderately 


deep  gratings,  g/d 


where  it  has  been  observed  for  the  latter 


that  as  g/d  increases  new  SP  dispersion  branches  pop  out  in  pairs  at  the 

w/ojp  «  1  //2  line.  As  g/d  keeps  increasing,  one  branch  of  the  pair  is 

pushed  upward  and  the  other  branch  is  pushed  downward  to  be  followed  by 

the  emergence  of  a  new  pair  of  branches.  Since  the  available  frequency 

interval  is  confined  both  from  below  (w/u>  -  0)  and  from  above  (w/w  - 

P  P 

1),  the  accumulation  of  the  branches  into  bands  as  6  -*■  ®  may  be 

expected.  However,  the  particular  band  gap  and  widths  we  derive  [Eqs. 

(4.4)  and  (4.5)]  may  be,  in  part,  an  artifact  of  our  approximations.  SP 

dispersion  bands  have  also  been  derived  for  an  infinitely  large  stack  of 

9 

thin  metallic  films. 


aL*-  -A.  *.  A* »  1-  •  -■  -  1*^'  - 


We  also  evaluate  a  simple  f igure-of-merit  representing  the  local 


field  enhancement  distribution.  The  local  field  enhancement  R  is 
defined  as  R  -  I  E/E .  |2,  where  E  is  the  total  electric  field  and  E. 

•  1  n  n  i  ’  1 


is  the  incident  electric  field. 


The  value  and  distribution  of  R  are 


of  prime  importance  in  determining  optical  and  chemical  processes  near 
the  grating.11  To  remain  within  reach  of  a  simple  analytical  result,  we 
evaluate  only  R(T/B)  -  |R(T)/R(B)|  ,  where  the  points  T  (top)  and  B 
(bottom)  are  defined  in  Fig.  1,  and  the  incident  light  is  normal  to  the 
grating.  The  result  [Eq.  (5.3)]  indicates  the  existence  of  two  distinct 
regimes,  i.e.,  X  >>  d  and  X  <<  d,  where  X  is  the  wavelength  of  the 
incident  light.  The  R(T/B)  ratio  in  the  former  is  smaller  than  in  the 
latter,  and  in  both  regimes  R(T/B)  -  1/&.  These  findings  are  understood 
in  terms  of  the  following  physical  picture:  The  incident  field  induces 
the  electrons  in  the  metal  into  an  oscillating  motion,  and  the  electrons 
respond.  Since  the  grating  is  deep,  the  motion  of  the  electrons  at  the 
grating's  peak  is  inhibited  by  the  grating's  boundaries,  and  much  less 
so  at  the  troughs'  bottoms,  and  hence  R(T/B)  <<  1.  The  distinction 
between  the  X  >>  d  and  X  <<  d  regimes  can  also  be  understood  in  the  same 
vein:  The  response  of  the  electrons  at  the  bottoms  is  inhibited  by  the 
grating's  boundaries  and  hence  depends  only  mildly  on  X.  The  response 
of  the  electrons  at  the  peak's  top,  however,  is  strongly  dependent  on  X: 
For  X  >>  d,  these  electrons  respond  very  poorly  due  to  the  confinement 
of  the  grating's  boundaries,  while,  as  X  ever  decreases,  the  top 
electrons  can  respond  ever  more  efficiently.  Consequently,  the  ratio 
R(T/B)  is  larger  for  X  <<  d  than  for  X  >>  d. 
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The  paper  is  organized  as  follows:  We  introduce  in  Section  2  the 

3 

dressed  Rayleigh  expansion  in  the  context  of  the  SG.  In  Section  3  we 
present  the  central  approximations  which  enable  a  simple  solution  for 
the  light  scattering  problem  in  the  asymptotic  limit  &  -*■  *.  The  results 
of  Section  3  are  then  applied  in  Section  to  analyze  the  SP  dispersion 
relations  and  in  Section  5  to  evaluate  the  local  field  enhancement 
distribution  for  a  particularly  simple  configuration.  Section  6  is  a 
discussion  and  summary. 


2.  Dressed  Rayleigh  Expansion  and  Sinusoidal  Grating 

3 

We  introduce  in  this  section  the  dressed  Rayleigh  expansion  and 
the  pertinent  formulas  of  the  SG.  As  has  been  demonstrated  elsewhere, 
the  dressed  Rayleigh  expansion  converges  for  an  SG  arbitrary  (5  and 
therefore  provides  a  convenient  framework  to  analyze  the  g  ■*  *  limit. 

All  expressions  in  this  section  are  exact.  The  notations  throughout  are 
defined  in  Fig.  1 

3 

The  dressed  Rayleigh  expansion  for  the  electromagnetic  fields  is 


En(x,z)  -  l  {Tn(£)p.  _U) 
u  l—  u’ 


i[kgX  -  Wq(£)(z  +  g)] 


+  a0U)  P0,*U)  e 


i[kgX  +  W1  U)  (z  +  g)  ] 


E^x.z)  -  l  Y2U)  p1  _(£)  e 


itk^x  -  W1(£)(z  -  g)] 


(2.1. a) 
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where 


Sq(x.z) 


.ffn  -  ”  itk  x  -  W  U)(g  ♦  z)] 

^  b  1  lX.U>  e  4  0 

i[k.x  +  Wn(fc)(g  ♦  z)] 

*  oQU)  e  1  } 


B1  (x,z) 


k  a  J  V*> 


itk^x  -  1W1 ( £.)  (g  -  z)] 


YnU)  -  60  n  En  e 

0  Ai  |  0  0  |  P 


♦iW0U)g 


,  (2.1 . b ) 


(2.2) 


k(i)  *  /e,  k,  k  «  u^/c,  k_  =  2n/d 
1  U  vj 


2  2  1/2 

WaU)  »  [k  (i)  -  ,  k^  -  k((  ♦  £kQ 


pi  ±u)  *  ^k£z  +  w1(£’)x]*  s  *=  X  x  z 


for  i  -  0,1 


(2.3) 


In  (2.1)  -  (2.3)  the  subscripts  of  the  fields  refer  to  the  domains  (0  or  1), 

-io)0t 

the  time  dependence  of  the  fields  is  assumed  to  be  e  ,  k((  is  the 

A  A  A 

■+ 

x-moraentum  component  label  of  the  field,  x,  z  and  s  are  unit  operators 
in  the  obvious  three  directions,  is  the  arbitrary  amplitude  of  the 
p-polarized  incident  light  normal  to  the  grating's  grooves,  and  we 
always  choose  Im[W^(l)]  or  Re[W^(£)]  £  0.  The  height  g  (positive)  is 
half  the  maximum  z-extension  of  the  grating  (Fig.  1),  and  the  central 


parameter  6  is  again 


e  *  gkQ  . 


(2.4) 


The  dressed  Rayleigh  expansion  coefficients  satisfy 


1  Vt  Vl)  ■  u(l") 


l  N  y AD  -  v(m)  , 

l*-m  m 'l  1 


(2. 5. a) 
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where  the  matrix  M,  for  instance,  can  be  dubbed  as  the  "Fresnel  matrix". 
It  is  given  by 


i ( W  Cm)  +  VL(l)g  _ 

M  .  “  e  M“  ,  u(m) 

m  v  x.  e  » Xr 


iW1(m)g  B 
e  p  (m) 


i(W  (m)  +  W  U))g  R 
N  -  e  N  n  »  v(m)  -  e 

HI  t  X.  ffi  |  X. 


lW0On)g  B 

v  (m)  , 


(2.5. b) 


and 


W  (£)W  (m)  +  kk 

M  0  -  -  u-rn - rj-r— r . -  i”"*  J  #[g(Wft(l)  -  W.(m))] 

m.t  Wq(£.)  ~  (m)  m- x.  0  1 

W  (t)W. Am)  *  k  k  - 

Nm,«.  “  W^JL)  -  W0(m)  1  Jm-£l-g(  _W1 U)  +  Wo(m))] 

-W  (0)W  (m)  ♦  kk 

^  (ra)  *  -"^(0)Vw^)-J  1  V-g(W0(0)  +  Vm))]E0p 


2  e_  e . 

V  (m)  -  ■■■■  ■  ~  Wn(m)6r_  _  E. 


c1  ~  e0  0  ni,0  ^0p  * 


(2.6) 


where  J  denotes  the  Bessel  function  of  order  v. 
v 

3 

It  has  been  shown  elsewhere  that  the  set  of  equations  (2.5)  can  be 
safely  truncated  at  a  maximum  | l | -value  L  of  the  order  5,  since  the 
matrix  elements  of  M  decay  exponentially  with  l  or  id  and  become 
insignificant  beyond  that  point.  A  more  detailed  examination  reveals 
that  the  important  matrix  elements  of  M  are  distributed  in  two  lobes: 
one  lobe  is  for  m  'small'  and  |£|  'large'  and  the  other  for  i.  'small' 
and  |m|  'large'.  The  diagonal  terms  diminish  exponentially  with  2.  »  m 
faster  than  either  lobe's  length. 

Once  the  set  of  equations  (2.5)  is  solved,  all  properties  of  the 
scattered  fields  are  known.  The  difficulty  when  B  ®,  however,  is  that 
the  relevant  matrix  M  (or  N)  is  very  large  and  highly  non-diagonal.  We 
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introduce  in  the  next  section  an  asymptotic  approximate  solution  of 
(2.5)  and  then  explore  its  implications  with  regard  to  the  SP  dispersion 
relation  and  light  scattering. 

3.  Approximate  Solution  for  6  >>>  1 

To  discuss  the  fields  scattered  by  the  grating,  it  is  necessary  to 
invert  the  Fresnel  matrix  M  [Eq.  (2.5)],  which  has  an  effective 
dimensionality  on  the  order  of  6  when  6  >>  1 .  Obviously  to  gain  insight 
we  need  to  introduce  further  simplifying  approximations.  The 
construction  of  a  simple  approximation  to  M  1  is  the  content  of  this 
section.  The  sequence  of  approximations  below  that  lead  to  the  final 
result  [Eqs.  (3-8)-(3. 10)3  assumes  that  B  >>>  1,  i.e. ,  very  deep 
gratings.  At  some  points  we  resort  to  additional  approximations 
motivated  by  the  gained  mathematical  simplicity.  We  believe  that  the 
plausibility  of  the  results  derived  in  Sections  4  and  5  based  on  the 
approximate  M  1  vindicate  the  approximations.  To  gain  perspective  on 
the  B  >>  1  regime  considered  here,  we  give  the  analysis  pertaining  to 
shallow  gratings  in  Appendix  A. 

An  attractive  feature  of  M,  or  M  ,  is  its  partial  separability  in 
the  "m"  and  "l"  indices  [Eq.  (2.6)].  As  will  become  obvious  shortly, 
separable  matrices  are  convenient  mathematically.  We  are  therefore 
motivated  to  devise  a  separable  approximation  to  the  other  factors  in 
(2.6),  in  particular  to  the  Bessel  function  factor.  This  turns  out  to 
be  possible  provided  |m|,  |S.|  and  | m—  JL )  are  very  large  compared  to 
unity,  which  is  the  case  for  most  of  the  elements  of  the  highly  non¬ 


diagonal  matrix  M. 
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n(B)  =  /l  +  B2  +  ln[B/(]  +  /I  +  B2)]  .  (3.l.b) 

The  above  is  valid  for  |m|  ,  |  £. [ .  ( m-£ f  >>  1.  It  can  be  easily  checked  by 

comparison  with  (2.6)  that  (3.1)  provides  a  surprisingly  good 

approximation  even  for  the  diagonal  element  . 

(S) 

Note  that  M  has  a  separable  form  except  for  the  denominator 

3/2 

[ Wq ( £ ) - ( m) ]  .  This  is  to  be  expected  since  a  fully  separable  form 

(S) 

for  M  does  not  possess  an  inverse.  To  maintain  the  simplicity  of  a 

separable  approximation  on  the  one  hand  and  yet  to  secure  the  existence 

of  a  simple  inverse  on  the  other  hand,  we  introduce  the  following 

additional  approximation:  Consider  first  the  off-diagonal  matrix 
(S) 

elements  of  M  .  These  elements  are  grouped  into  two  lobes  (Section 

2),  i.e.,  when  |m|  is  'small'  and  |£|  is  'large',  and  vice  versa.  When 

|m|  is  'small'  and  |£|  is  'large',  for  example,  the  fastest  |£|-changing 

( S) 

factor  in  M  is  the  exponential  [Eq.  (3.1)1.  Consequently,  the  above 
mentioned  now  separable  denominator,  which  changes  with  j  £ |  much  more 


!VVW 


TTTtTv 
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slowly,  acts  approximately  as  a  normalizing  factor.  It  can  be  therefore 

approximated  by  a  suitably  chosen  constant,  denoted  by  C1 .  For  the  sake 

■3/2 

of  simplicity,  we  further  approximate  [WQ(  (m)]'  -  for  both  the 

/  c  \ 

|m|  <<  |i|  and  |fc|  <<  |m|  elements.  The  diagonal  elements  of  Mv 

comprise  the  complementary  (minority)  group  of  elements  to  the  off- 

diagonal  elements  just  discussed.  In  this  case,  the  above  argument 

(S) 

concerning  the  fastest  | £ | - | m| -changing  factor  in  M  is  still  valid. 
However,  the  corresponding  constant  approximating  the  non-separable 
denominator  is  obviously  quite  different,  i.e.,  [WQ(  O-V^  ( 03  «  C 2» 

The  foregoing  discussion  suggests  therefore  the  approximation 


M  ~  M^S5 
m,  l  m,  S. 
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(3.3) 


and  C1  has  been  defined  in  the  discussion  above. 

The  diagonal  elements  in  Eq.  (3.2)  depend,  of  course,  on  the 

choice  of  C.and  C  .  However,  regardless  of  the  particular  choice,  due 

ig[W0U)  ♦  W1  (m)  ] 

to  the  convergence  factors  e  ,  the  number  of 

significant  diagonal  elements  is  of  the  order  of  unity.  Therefore, 
there  are  only  very  few  significant  diagonal  terms  in  comparison  with 
the  number  of  off-diagonal  terms,  yet  if  they  are  neglected,  Eq.  (3-2) 
ceases  to  have  an  inverse.  We  have  therefore  a  situation  where  a  few 
diagonal  terms  are  not  important  in  terms  of  physical  results,  yet  due 
to  the  particular  mathematical  structure  (separability)  it  is  necessary 
to  keep  them  to  assure  the  existence  of  an  inverse.  Under  such 


conditions,  it  is  probably  immaterial  as  to  how  the  diagonal  terms  are 


chosen,  provided  the  calculated  observable  does  not  depend  on  the 

specifics  of  the  choice.  We  therefore  simplify  Eq.  (3.2)  further  by  the 

«.  B 

approximation  -  D,  where  D  is  a  very  small  number  of  the  order  e  or 
less.  Again,  to  be  consistent  with  the  arguments  above,  we  can  use  the 
ensuing  approximation  only  to  calculate  observables  that  do  not  depend 
on  D,  such  as  those  analyzed  in  the  next  two  sections.  We  would  like  to 
emphasize  that  the  -  D  approximation  is  not  essential,  i.e.,  it  is 
possible  to  invert  Eq.  (3-2)  as  it  stands.  The  substantial 
simplification  gained  by  adopting  this  approximation,  however,  warrants 
an  examination.  Therefore,  Eq.  (3.2)  is  replaced  by 


M  .  -  D<5  +  V  . 

m,  8.  m,  1  m,  l 


(3.4) 


-3 


where  D  is  a  very  small  number  of  the  order  e  . 

The  inversion  of  Eq.  (3-4)  can  be  obtained  now  in  a  straightforward 
manner.  Note  from  Eq.  (3. 3)  that  V  is  comprised  of  four  separable 
terms, 


V  .  -  l  a  (m)  b.U)  , 
">U  ^  J  j 


(3.5) 


where  the  factors  a^  and  b^  are  extracted  from  Eqs.  ( 3 * 1 • b )  and  (3-3) 
[see  Appendix  C,  Eq.  (C.2)].  The  form  of  Eq.  (3.5)  has  the  attractive 
property  that  V  to  the  n-th  power  is  given  by 


(Vn) 


m,  8, 


1  a  (m)  [An_1  ]  .  b.U)  , 

i, j-1  J  J 


(3.6) 


where  the  auxiliary  interaction  matrix  A  is  of  dimension  4  with  the 
elements 


A,  -  l  b,U)a4U)  . 


(3.7) 


The  auxiliary  matrix  A  (Appendix  C)  plays  a  central  role  in  all 
subsequent  manipulations.  Consequently,  expanding  the  inverse  of 
Eq.  (3.*0  in  a  geometric  series  and  using  Eq.  (3*6)  yields 


D 


^  l  a  (m)([1  +  4a]  ])  b  (£) 
D  i,j-1  D  i.J  J 


1_ 

D 


l  a  (m) [A  1  ]  b  U)] 

i,j-1  1  1,J  J 


(3.8) 


The  second  line  in  Eq.  (3.8)  follows  since  <|a|>/D,  where  <  J  A ] >  is  a 
typical  average  value  of  the  elements  of  A  [see  ( C . 18)  and  (C.20)]  and  D 
is  an  extremely  small  number. 

The  approximate  Inverse  (3-8)  is  a  key  result  of  this  work.  The 
unknown  parameter  D  appears  as  a  multiplicative  factor,  implying  that 
(3.8)  can  be  used  only  to  evaluate  ratios  so  that  D  drops  out  from  the 
final  result.  Note  that  (3.8)  embodies  the  reduction  of  an  inversion 
problem  of  a  B  *  B  matrix  to  the  inversion  of  a  ^  ^  matrix. 

Based  on  the  analytical  expressions  for  the  matrix  elements  of  A 
(Appendix  C) ,  the  inverse  A  1  can  be  further  simplified.  The  key  remark 
is  that  when  B  •»  ®,  the  matrix  elements  denoted  by  "a"  in  Table  1  are 
substantially  smaller  than  the  rest  of  A  [see  (C.11)  and  following 
discussion]  and  therefore  can  be  neglected.  Consequently, 

det(A)  =  det(Av)  det(A.) 

—  —  T  —0 


det(AY)  -  A^  2  A^  ^  -  A1  ^  A^  2 

det(A^)  =  y  A^  ^  ~  A^  3  A^  ^  , 

(3.9 

where  the  symbol  "det"  indicates  a  determinant  ,  and 
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+  mixed  y  and  6  terms 

The  simplicity  of  (3.9)  and  (3.10)  in  conjunction  with  (3.8)  is  most 
gratifying.  It  indicates  an  expectation  that  in  the  limit  B  >»  1,  as 
in  the  B  <<  1  limit,  the  physics  should  simplify.  All  that  is  needed  to 
obtain  M  1  are  eight  matrix  elements,  combined  in  the  simple  form  of 
(3.10)! 

4.  Surface  Plasmon  Dispersion  Relation 

Surface  plasmon  excitations  are  collective,  surface-attached 

oscillations  of  the  metal's  electron  gas,  well  studied  in  the  context  of 

7  S 

flat  metallic  surfaces  and  shallow  gratings.  ’  For  the  case  of  a 

shallow  grating,  the  dispersion  relation  is  comprised  of  many  branches 

(for  a  lossless  metal,  as  assumed  here  throughout)  which  are  arranged  in 

pairs  above  and  below  the  line  k/k  -  /2,  where  k  -  u  /c  and  u  is  the 

P  P  P  P 

volume  plasmon  frequency.  In  this  section  we  derive  the  dispersion 

relation  for  the  B  ■*  •  limit,  i.e.,  pertaining  to  very  deep  gratings. 

The  emergence  of  pairs  of  dispersion  relation  branches  when  B  **  0 

can  be  understood  using  arguments  employed  in  interpreting  the  opening 

1  2 

of  a  band  gap  for  electrons  moving  in  a  periodic  ionic  potential.  In 


u 


the  latter  ease  there  are  two  Interfering  combinations  of  unperturbed 

electron  wave  and  the  umklapp-reflected  wave.  Since  one  combination  is 

centered  at  the  ions'  sites  and  the  other  between  ions,  the  associated 

energy  of  these  two  combinations  are  split  and  a  band  gap  opens.  In  the 

present  case  the  Bragg  reflections  play  the  role  of  the  umklapp 

processes  and  the  unperturbed  electron  waves  —  by  the  flat-surface  SP 

mode  (see  Appendix  A)and  hence  the  pairs  of  dispersion  branches.  We 

expect  this  type  of  dichotomy  to  persist  for  &  ■*  ®.  The  remaining  issue 

is  therefore  how  do  these  pairs  of  branches  accumulate  for  &  -*•  ®,  given 

that  the  available  frequency  domain  is  bounded  from  below  by  u)/u>p  -  0 

and  from  above  by  ui/w  »  1 . 

P 

By  definition,  the  SP  excitations  decay  exponentially  in  both  z- 
directions.  Fig.  1,  i.e.,  these  are  excitations  for  which  E^  -  0  in 
(2.6).  Consequently,  to  have  a  non-trivial  solution  of  the  Fresnel 
equations  (2.5),  the  condition 

det(M)  -  0  (4.1) 

determines  a  function  of  k(0)  versus  k ,,  provided  that  -k_/2  £  k,.  £  k_/2, 

II  U  II  u 

7 

i.e.,  k||  is  confined  to  the  first  Brillouin  zone.  Condition  (^.1)  is 
tantamount  to  finding  the  points  where  M  1  does  not  exist  (poles). 

Hence,  in  the  context  of  the  approximate  inverse  (3.8),  Eq.  (H.1) 
transcribes  into  det(A)  -  0.  The  determinant  of  A,  however,  is  given  by 
(3.9)  in  the  leading  asymptotic  order.  Consequently,  since  det(A.)  *  0 
outside  the  light  cone  [see  Eq.  (C.22)],  the  SP  dispersion  relation  is 


det(A^) 


A 


1,2 


A1 . 4  A3.2 


0 


(4.2) 


-> 


simply 


Inserting  into  (4.2)  the  explicit  expressions  for  the  matrix  elements 
( C. 20)  yields 

[3k2(0)  +  k2(1))(k2(0)  ♦  3k2(l)3 
3[ k2(0 )  +  k2(1)]2 

sin[L(a  +  Y)]sln[-  -]  -  sinL[(a  -  Y)3sin[a  ^3  5 

“  [ _ t _ £ _ y  (4. 3. a) 

sin[L(a  +  Y) 3sin[— -3  +  sin[L(a  -  Y)3sln[-a-^— -3 

where  k(i)  «  k  from  (2.2)  and 

Y  -  2iTk  /k 

U 

a  -  ^-(-[k2(0)  ♦  k2(1)3)1/2  .  ( -4  - 3 - 1> ) 

KG 

In  (4. 3. a)  L  is  a  very  large  number  of  the  order  &  [see  the  discussion 

after  Eq.  (C.12)3,  the  precise  value  of  which  is  immaterial  (see  below). 

The  analysis  of  (4.3)  is  quite  simple.  Consider  first  the  k (0 )  £ 

1  /2 

kp/[1  +  Cq3  regime,  where  a  and  obviously  Y  of  (4.3.b)  are  real  and 

hence  the  RHS  of  (4.3. a)  is  real  and  positive.  It  can  be  easily  shown 

that  for  (4.3. a)  to  have  a  solution  it  is  necessary  that 

k  <  k  /[3eA  +  1 31 /2.  When  k(0)  2  k  /[I  +  Gn31/2  the  a-parameter  is 
p  u  p  u 

imaginary  and  hence  the  RHS  of  (4. 3. a)  is  real  and  negative.  Simple 
considerations  show  again  that 

to  make  the  LHS  of  (4.3. a)  negative  it  is  necessary  that  k  <  [•=■— - 3^ 

5  +  e0 

in  order  for  a  solution  to  (4.3)  to  exist,  it  is  necessary  that 
(assuming  eQ  -  1  for  simplicity) 

k  k  rr 

0SkS;r£,  ^  k  <  ~  k  .  (4.4] 

2  /2  2  p 

Condition  (4.4),  however,  is  also  sufficient.  Since  L  -  B  >>>  1,  the 
RHS  exhibits  oscillations  which  are  ever  increasing  in  amplitude  and 


'  ‘  ■ 
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MV 


& 
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frequency  as  $  is  increased.  Therefore,  any  k-value  within  the 
admissible  ranges  (4.4)  generates  a  solution  of  (4.3),  i.e.,  the 
dispersion  relation  constitutes  two  bands  (Fig.  2).  The  frequency  gap 
Aw  between  the  bands  is  according  to  (4.4) 


Aw  /2  -  1 
w_  2 


(*1.5) 


p-  > 


The  results  ( *4 . 4 )  and  (4.5)  describe  the  accumulation  of  dispersion 
branch  pairs  mentioned  in  the  beginning  of  this  section.  The 
interpretation  of  the  two  bands  is  identical  to  that  of  the  two  members 
of  a  dispersion  branch  pair:  The  two  bands  are  comprised,  respectively, 
of  modes  in  which  the  field  is  mostly  localized  either  in  the  troughs  or 
in  the  peaks  of  the  grating.  We  cannot  be  sure  at  this  point  to  what 
degree  the  particular  band  gap  (*1.5)  and  band  widths  are  an  artifact  of 
our  approximations.  At  the  very  least,  C4 . 4 )  and  (4.5)  indicate  domains 
of  high  density  of  dispersion-relation  curves. 

Note  that  the  dispersion  relation  (4.3)  or  (4.4)  does  not  depend  on 

g,  k„  or  k„.  The  independence  of  g  or  k„  can  be  understood  as  follows: 

Ci  l!  u 

If  6  -*•  ®  is  Interpreted  as  g  ♦  ■  and  d  finite,  then  obviously  the  result 
of  any  observable  cannot  depend  on  g  in  the  asymptotic  region.  By  the 
same  token,  the  limit  B  ■*  *  can  be  construed  as  keeping  g  finite  but 
letting  d  +  0,  which  leads  again  to  the  conclusion  that  the  result 
cannot  depend,  to  the  leading  order,  on  k*.  The  independence  of  the 

U 

band  structure  of  <  may  be  an  artifact  of  our  approximations.  The 
behavior  of  the  matrix  elements  of  M  near  the  light  cone,  for  instance, 
is  controlled  by  the  W^(fc)  sectors  [Eq.  (2.2)],  with  the  square-root 
feature  playing  an  essential  role  (Appendix  A).  On  the  other  hand,  in 


all  manipulations  employed  here,  we  expand  the  square  root  W^(£.)  as 
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*0 ' 


•■V 


V£)  =  iN“il  -  )L^Q-  • 


(H.6) 


which  destroys  the  regardation  effects  near  the  light  cone. 

5.  Local  Field  Enhancement  Distribution  for  Normal  Incident  Light 

As  a  second  application  of  our  central  approximation  (3-8),  we  turn 
now  to  the  calculation  of  the  local  field  enhancement  distribution  at 

the  grating's  surface.  The  local  field  enhancement  R(x,z)  is  defined  as 

1  2 


R(x,z) 


'0(x  ,z) 


It  measures  the  collective  response  of  the 


"Op 


electrons  in  the  metal  to  an  incident  electromagnetic  field,  and  is  of 
obvious  importance  in  determining  optical  and  scattering  processes  near 
the  grating.11  To  simplify  the  analysis,  we  pick  the  simplest 
configuration.  We  assume  normal  incidence  of  the  p-polarized  light  and 
consider  the  ratio  of  the  local  field  enhancements  at  the  top  of  the 
peak  to  that  at  the  bottom  of  the  trough  (points  T  and  B  in  the  notation 
of  Fig.  1).  This  ratio,  denoted  by  R(T/B) ,  is  a  convenient  figure-of- 
merit  describing  the  local  field  enhancement  distribution.  Since  R(T/B) 
is  a  ratio  of  two  fields,  the  unknown  constant  D  in  (3.8)  does  not  enter 
the  final  expressions. 

Combining  (3.8),  (3-10),  (2.5)  and  (2.1),  it  is  straightforward  to 
obtain  for  the  electric  field  in  region  0  (Fig.  1) 

i[knx  -  W_(0 )z]  - 


£q(x,z)  »  YQ(0)e 


po,- 


(0) 


.  "  -  i[k.x  ♦  WnU)z] 

♦  ^[{  I  PQ  +U)  e  mU) 

t»-00  ’ 


-II  I  P0  +U)  e 
j ,  n-1  ’ 


i[k  x  +  W.U)z] 


,U)}(  A-1).  n{  l  b  (m)  y(m) }  ] 
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The  decomposition  of  EQ(x,z)  in  the  last  line  of  (5.1)  reflects  the 
structure  of  (3.10):  The  first  term  in  (3-10),  involving  the  matrix 
elements  A1  2>  etc.  (or  class-Y  elements  in  the  terminology  of 

Appendix  C  and  Table  1)  gives  rise  to  field  components  that  vary  very 
rapidly  with  small  changes  in  either  k,  k  or  k„,  etc.  This  feature  is 

ll 

manifested  in  the  explicit  expressions  for  these  matrix  elements  [Eq. 

(C.20)]  and  has  been  exploited  already  in  conjunction  with  the  SP 

I  dispersion  relation  (Section  4).  The  corresponding  component  in  the 

total  field  is  denoted  by  Ef(x,z),  where  the  subscript  "f"  indicates 

"fast".  The  second  term  in  (3-10)  involving  the  elements  A2  l ’  A2  3* 

|  etc.  (class-6  in  the  terminology  of  Appendix  C  and  Table  1)  changes 

slowly  with  small  changes  in  k,  k^ ,  etc.  and  gives  rise  to  the  "slow" 

component  In  the  total  field,  denoted  by  E  (x,z)  in  (5.1).  The  incident 

s 

|  field,  denoted  by  Einc(x,z),  is  the  first  term  on  the  LHS  of  (5.1).  The 

decomposition  (5.1)  of  the  total  field  into  a  fast  oscillating  component 
£  (x,z)  and  a  slow  envelope  E  (x,z)  is  expected  due  to  the  many 
I  interfering  components  (Bragg  reflections)  of  which  EQ(x,z)  is 

comprised. 

In  terms  of  the  decomposition  (5.1),  we  identify  the  ratio  of  the 
local  field  enhancements  at  points  T  and  B  in  Fig.  1  with  the  field 
envelope,  or  average: 


R(T/B) 


i 


E.  (T)  ♦  E  (T) | 
inc  s 

2 

«s(T) 

E .  (B)  ♦  E  (B) 

Inc  s 

es(B) 

(5.2) 


where  the  arguments  "T"  and  "B"  in  (5.2)  indicate  the  coordinates  of 
points  T  and  B.  Straightforward  yet  somewhat  lengthy  algebra 
(Appendix  D)  gives 
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R(T/B)  *  ^  for  r  «  1 


and  ( eQ  *  1  for  simplicity) 


(5.3) 


(5.10 


The  complete  expressions  are  given  in  Appendix  D,  whereas  (5.3) 
represents  an  ordei — of-magnitude  estimate. 

The  result  (5.3)  is  interpreted  in  the  following  way:  The  fact 
that  there  are  two  regimes,  crudely  defined  by  X  <<  d  and  X  >>  d,  can  be 
understood  in  terms  of  the  electrons'  relative  response  to  an  incident 
field  at  the  top  of  the  peaks  and  the  bottom  of  the  troughs.  When 
X  >>  d  (r  <<  1),  the  electrons  at  the  peaks'  tops  respond  poorly  by 
comparison  to  the  electrons  at  the  bottom  of  the  troughs,  simply  because 
the  former  are  confined  to  the  metallic  fingers  of  the  grating.  On  the 
other  hand,  when  X  <<  d  (r  >>  1),  it  is  obvious  that  the  response  of  the 
electrons  at  the  peaks'  top  is  impeded  to  a  much  lesser  degree  by  the 
boundaries  of  the  grating's  fingers,  while  the  electrons  at  the  troughs’ 
bottoms  respond  in  about  the  same  manner  as  in  the  X  >>  d  regime.  This 
consideration  implies  that  R(T/B)  is  smaller  for  X  >>  d  than  for  X  <<  d, 
which  is  born  out  by  (3-3).  The  fact  that  in  both  regimes  R(T/B)  -  6  ^ 
<<  1  is  in  keeping  with  the  above  picture:  Since  the  grating  is  deep, 
for  normal  incidence  it  is  always  easier  for  the  electrons  to  respond 
when  they  are  located  at  the  bottoms  of  the  troughs  than  when  they  are 


at  the  tops  of  the  peaks.  Note  that  these  considerations  do  not  depend 

on  g,  since  they  pertain  to  relative  magnitude.  It  is  interesting  to 

notice  that  the  same  trends  are  present  in  numerical  studies  of  shallow- 
10 

gratings. 

6.  Discussion  and  Summary 

The  main  message  of  this  work  is  that  for  the  very  deep  sinusoidal 
grating  limit  (6  >>>  1),  the  physics  simplifies  considerably:  By 
introducing  a  sequence  of  reasonable  approximations,  the  problem  of 
inverting  the  Fresnel  matrix  (2.6)  of  order  &  has  been  reduced  to  the 
calculation  of  eight  matrix  elements  (3.10).  The  simplicity  of  the 
outcome  is  in  keeping  with  a  related  study  of  the  surface-plasmon 

9 

dispersion  for  an  infinite  stack  of  metallic  thin  films. 

With  regard  to  the  surface-plasmon  dispersion  relation,  we  find  two 

bands,  separated  by  a  frequency  band  gap  given  by  Aw/wp  »  0.2.  These  two 

bands  represent  two  classes  of  localized  SP  excitations:  Those  for 

which  the  fields  are  maximized  in  the  troughs  and  those  where  the  fields 

are  maximized  inside  the  peaks.  The  latter  class  of  modes  correspond  to 

the  upper  band  since  they  represent  an  enhanced  oscillation  of  the 

electrons.  This  interpretation  is  in  keeping  with  the  interpretation  of 

dispersion  branch  pairs  already  observed  for  shallow  gratings.  It  is 

also  consistent  with  the  analogy  of  the  two-adjacent-band  structure  of 

1 2 

electrons  moving  in  a  periodic  potential.  With  regard  to  the  local 
field  enhancement  distribution,  we  find  that  the  field  at  the  bottom  of 
the  troughs  is  always  larger  than  at  the  tops  of  the  peaks.  However, 
the  relative  magnitude  of  the  "top/bottom"  enhancement  exhibits 
different  behavior  in  the  A  >>  d  and  A  <<  d  regimes.  The  existence  of 


these  regimes  reflect  primarily  the  response  of  the  electrons  at  the 
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peaks'  tops.  The  electrons  there  are  restricted  In  their  motion  by  the 
grating's  boundary,  and  therefore  respond  poorly  when  X  >>  d  and  much 
better  when  X  <<  d. 

The  simple  outcome  of  this  study  will  hopefully  stimulate  further 

studies  in  the  same  vein.  To  clarify  to  what  degree  the  specifics  of 

our  results  do  not  depend  on  an  artifact  of  the  approximations,  it  would 

be  desirable  to  attempt  some  refinements,  e.g.,  the  straight  inversion 

of  the  quasi-separable  approximation  (3-1)  or  (3-2).  Another 

possibility  is  to  carry  out  large-scale  exact  calculations  of  the  SG 

3 

using  the  dressed  Rayleigh  expansion  as  described  in  Section  2.  In 
particular,  it  is  interesting  to  substantiate  the  quantitative  values  we 
obtain  for  the  width  and  gap  of  the  SP  bands  [Eq.  (4.4)  and  (4.5)]. 

These  results  hopefully  will  generate  interest  also  for  corresponding 
experimental  work. 

Acknowledgments 

We  are  grateful  to  A.  C.  Beri  for  his  assistance  in  preparing  the 
figures  and  to  J.  Rouse,  R.  Ferron  and  L.  W.  Edwards  for  their  typing  of 
the  manuscript.  This  work  has  been  supported  by  the  Office  of  Naval 
Research  and  the  Air  Force  Office  of  Scientific  Research  (AFSC),  United  States  Aii 
Force,  under  Contract  F49620-86-C -009 .  The  United  States  Government  is 
authorized  to  distribute  reprints  notwithstanding  any  copyright  notation 
thereon. 


-  .  ■  ■  «•  v  - 
•  .  •  „  •»  ..  »  *  *  «  ^  u  ■'  s  *■  +  “*«  •  ‘  ■* 


Appendix  A:  Plasmon  Dispersion  Relation  for  a  Shallow  Grating  Surface 

To  highlight  the  particular  features  of  the  B  <<  1  regime  and  to 

appreciate  the  role  of  the  various  factors  in  (2.6),  we  analyze  here  the 

shallow  grating  (B  <<  1)  SP  dispersion  relation.  The  results  are  not 
5  7  8 

new,  '  ’  yet  they  are  obtained  here  in  a  direct  and  transparent 
manner. 

Consider  first  the  case  of  a  flat  surface,  i.e.,  g  -  0  in  the 

1  4 

notation  of  Fig.  1.  Since  J^(z=0)  «  6^  Eq.  (2.6)  yields  m -l  and  the 
matrix  M  is  diagonal.  Consequently,  for  the  determinant  of  M  to  vanish, 
Eq.  (4.1),  one  or  more  of  the  diagonal  elements  must  vanish.  Hence,  the 
dispersion  relation  reads 


WqUJWjU)  +  k^  -  0 

or,  by  inserting  the  definitions  (2.3),  (A.1)  yields 

.2  .2  e0e1 
k .  *  k  -  . 

1  Vci 


(A.1) 


(A. 2) 


For  a  flat  surface  d  -*  ®,  i.e.  k  »0,  so  that  k.  k  and  (A. 2)  is  the 

U  X. 


well-known  dispersion  relation. 


7,8 


Thus  the  numerator  in  (2.6) 


involving  Wrt(£)W,  (m)  +  k.k  entails  the  retardation  effects,  namely, 

0  1  x.  m 

that  the  SP  dispersion  relation  cannot  cross  the  boundaries  of  the  light 


cone. 


Consider  now  the  case  of  a  very  shallow  grating,  i.e.,  g  is  finite 

yet  6  -  gk  <<  1  .  Here  the  argument  of  the  Bessel  function  factor  is 
0 

1 4 

small,  yet  finite.  Therefore,  since 

lim  J  (z)  ;  (-)n  for  v*0 
n  n  v 

z+0 

r  1  for  v  =  0  (A. 3) 

to  first  order  in  8,  the  only  non-diagonal  terms  that  contribute  are 

when  |  m- 2.  |  <  1.  Thus,  M  is  tri-diagonal,  with  a  dominating  diagonal  and 


the  two  first-off-diagonal  of  order  6.  The  B  -  0  and  B  <<  1  examples 
discussed  so  far  clearly  indicate  that  the  Bessel  function  factor  in 
(2.6)  determines  the  amount  of  mixing  of  the  m  and  £  Bragg  reflections. 

The  non-vanishing  denominator  of  (2.6)  affects  only  the  normalization 
and  the  fall-off  in  m,£  rate  of  H 

m  i  *> 

To  determine  the  conditions  for  (4.1)  to  hold  when  B  <<  1  is  no 

simple  task  in  general.  However,  for  special  (k,k^)  combinations  the 

problem  can  be  easily  solved.  Let  us  thus  consider,  in  particular,  the 

boundary  of  the  first  Brillouin  zone  k ,,  »  k^/2  and  the  crossing  of  the  £ 

II  u 

-  0  and  £  »  -1  dispersion  curves  given  by  (A. 2).  If  there  were  not 

off-diagonal  terms  in  M,  the  crossing  of  the  curves  at  k  ,  =  k./2  would 

—  It  u 

define  a  (k,k^)  point  for  which  (4.1)  is  satisfied.  For  finite  g, 
however,  the  two  relevant  off-diagonal  elements  M  and  M  change 

U,  - J  ”1  |U 

the  eigenvalues  of  M.  The  condition  for  the  perturbed  eigenvalues  of  M 
to  vanish  [so  that  (4.1)  is  satisfied]  is  therefore  that  the 
corresponding  determinant  vanishes: 

»  0  (A. 4) 

k()  =  kQ/2  ,  B  <<  1 

Eq.  (A. 4)  is  an  equation  for  k  which  has  two  solutions,  whose 

difference  is  the  band  gap.  This  equation  can  be  easily  solved  provided 

that  at  the  Brillioun  zone  boundary  (k  =  k„/2)  we  are  allowed  to 

II  u 

neglect  retardations,  i.e.,  kQ  >>  k(0),|k(1)j.  In  this  case,  WQ(0 )  - 
W  ( -1 )  -  W  (0)  -  Wn(-1)  -  ik  /2,  and  some  simple  algebra  gives  the  gap 
frequency  Aw  as^ 


WVVV 


'  *v W  -  .'v. 


a-3 


I 


Appendix  B:  Separable  Approximation  for  the  Bessel  Function  Factor  In 

(2.6) 

We  derive  a  separable  approximation  to  the  Bessel  functin  factor  in 
the  matrix  M: 

-  !B-,) 
which  is  valid  for  8  >>  1  and  |i,|,|m|  >>  1.  The  second  line  in  (B.1) 
follows  from  (4.6).  We  first  analyze  (B.1)  by  considering  separately 
the  pertinent  asymptotic  domains  and  then  introduce  a  simplifying 
separable  expression. 

When  1  <<  |  Jt-m |  <<<  8,,  the  argument  of  the  Bessel  function  (B.1) 
is  much  larger  than  its  order.  Correspondingly,  the  leading  asymptotic 
term  is 

Vz)  =  cosU-tf+f)] 

V  -  m-i.,  z  “  i  B (  1 5L  |  - 1  m | )  r  g[wQ(  fi. ) -W1  (m)  ]  .  (B.2) 

When  8  -  |!*m|,  the  argument  of  the  Bessel  function  (B.1)  and  its  order 

are  large  and  proportional  each  other.  The  asymptotic  expression  in 

.  .  14 

this  regime  is 


1 


J  (vz)  z 
v  /2nv 


[1-z2] 


1/4 


v’i'(  z) 

e  + 


e(v-3/2) 


v(z) 


/I  -2 


+ 


ln[z/( 1 + 


(B.3) 


for  |arg  z  |  <  u  ,  v  ■+  » 

To  apply  (B.3)  to  (B.1),  we  consider  separately  the  three  possible  sub- 

domains:  (a)  2.,m  >  0  (b)  £,m  <  0,  (c)  |£,|  >>  1,  |m|  >>  1  and  m,i.  of 

opposite  signs.  Using1^  J  (z)  -  (-1)n  J  (z)  and  J  (-z)  «  (-1)nJ  (z),  a 
rr  n  n  n  n 

simple  a  , sis  gives 


1 


1 


A- 


[i ♦e2] 


i  a 


i-m| n( 6) 


n|  6) 


fcn[B/(1  + 


(BA) 


where  B+  refers  to  domain  (a)  and  B  to  domain  (b).  When  £.  and  m  have 
different  signs  [domain  (c)],  the  argument  of  the  Bessel  function  (B.1) 
is  much  smaller  than  its  order,  and  hence  the  corresponding  asymptotic 
expression  is1^ 


J  (z)  z  -  (f^-)  ,  v  «  m-i,  z  -  i&(  £  -  m  )  .  (B.5) 

v  rx —  2v  i  i  i  i 

/2ttv 

Consequently,  z/v  -  3  | ( £+m) /( t-m) j  <<  1  and  hence  B^  ^  *  0  for  most  of 
the  m,!t  combinations. 

Given  (B.2)~(B.5) ,  the  approximation 


=  /p 

V2i 


1 


1 


u.b2)'"1 


exp(i<{>)  +  exp(-iip)} 


*  *  4^  CWn(£)-W,(m)]  -  [^4]t 


(B.6) 


unifies  all  relevant  domains. 


Appendix  C:  Reduced  Auxiliary  Matrix  A 

We  calculate  some  and  estimate  the  remaining  matrix  elements  of  the 

auxiliary  matrix  A  defined  by 
00 

Ai  4  -  l  b  U)a  (i)  ,  i,j  -  1,  ..n  (c.  1 ) 

,J  J 

where,  from  (3-1)  and  (3.3), 

a1 (m)  -  W1 (m)  exp[i(W1(m)(g-n(B)/kG)-  J)] 

b2U)  =  cwQU)  exp[i(M0(l)(g*n(B)/kG))] 

a2(m)  *  W1  (m)  expfiCW^  (m)  (g+nU)/kG)  +  to  +  J)  ] 
b2(i)  *  CWQU)  exp[i(W0(£)(g-n(B)|kc)  -  *£)] 

a3(m)  =  km  exPt1(wi  (m)(g_n(8)/kG)  -  |j)] 
b3U)  -  Ck£  exp[i(W0U)(g+n(B)/kG))] 

a^(m)  =  km  exp[i(W1 (m) (g+n(B)/kQ)  +  to  +  J)  ] 

b^(£)  -  Ck£  exp[i(W0(i)(g-n(g)/kQ)  -  itl)  ]  (C.2) 

and  C  is  defined  in  (3-3).  Note  that  in  (C.1)  all  the  independence 
enters  via  the  combintion  k£  -  k(|  +  £kQ  [see  (C.2)].  Hence,  the 
summation  in  (C.1)  can  be  evaluated  by  using  the  Poisson  formula:15 

CO  CD  w 

l  G(y+U)  =  j  l  g(  2irn/ A)  e"2luny/A  g(t)  -  [  du  G(u)eitu  (C.3) 

i.*-®  m=-»  )  _oo 

with  the  identif ictaions  y+k^  ,  &+kG  and  u-+k£.  Furthermore,  since  for 

most  terms  in  (C.1)  |i|  >>  1,  it  is  justified  to  evaluate  A  .  using  the 

i » J 

asymptotic  expression  of  WjU),  Eq.  (4.6).  Keeping  the  leading  term  in 

(4.6)  is  not  sufficient  since  all  references  to  k(j)  enter  in  the  second 

term.  It  is  therefore  the  first  small  connection  to  A  ,  beyond  the 

*  *  J 

leading  term,  that  carries  all  the  physics. 


Inserting  (^.6)  into  (C.2),  it  is  easily  seen  that  all  Fourier 


transforms  g(u)  mandated  by  (C.3)  take  the  form 

f"  L*il  f“  i[f  (Kn)+lc,t]0 

giJ(t)  *  \_«kl  bi(k£)aj(k£)e  *  J«ct  We 

for:  i,j  -  1  .... 4 ;  n  -  and  n  -  a,B,T,6  .  (C.lJ) 

where 

F1(kl}  “  W0U)  W1U)  =  cl'ki  +  |[k2(0)*k2(  1 )  ] } 

F2(kjL)  -  C (~ik^|  k^|  ♦  |k2(0 ) sgn(k^)  { 

W  -  Ci-ikJkJ  ♦  |  k2(  l)sgn(k^) } 

F4^KjL)  “  C  {k2  }  (C.5) 

and 


fa  *  +  l  +  w0U)(g-n(6)/kc  ♦  W^Otg+nU)/^] 

J  +  2ig|kt|  -  igk2(1)-|^j 
f6  “  '  5  +  w0^)tg+h(S)/kG]  ♦  w1U)[g-n(e)/k(,] 

-  ~  +  2ig|'<i|  -igk2(0) 

ty  *  J  +  w0(£)[g+n(e)/kG]  ♦  w1(£)rg+n(e)3  ♦  u£ 

=  l  -  ^,/kQ  +  ♦  Mig|  k£j  -  ig[k2(0)  +  k2(  1 )  ]/|kJ 


f6  “  -  5  +  W0U)[g-n(e)/kG]  ♦  W1  ( £) Cg-n(  6)  |  kG3  -  nil 

♦  -  !l!i  ♦  J_  |k  |  .  1 [k2(0 )+k2(  1 )  ]  1 

k_  k.  Bk_  I  2.1  ilftu  TIT 


IT 

5 


T*7T  • 


(C.6 ) 


The  sixteen  Fourier  transforms  g  (t)  [Eq.(C.M)]  are  obtained  from  ai: 

1  »  J 

combinations  of  the  F^  and  the  f^.  The  explicit  correspondence  in 
(Table  1 )  is 


A-7 


1 

with 

CM 

CM 

< 

A2,4* 

A4 ,2  ’ 

A4 ,  4 

£  fe 

with 

Ai  ,1  • 

A1,3* 

A3.i* 

A3.3 

fY 

with 

A1 ,2* 

A1  ,4* 

A3.2' 

A3,4 

L  f  6 

with 

A2,1  * 

A2,3* 

A4,1  ’ 

A4,3 

(C.7 ) 


The  procedure  of  calculating  A  is  now  clear.  We  first  evaluate 

1  t  J 


the  transforms  (C.4)  for  all  t-values  mandated  by  the  Poisson  formula 
(C.3).  and  then  sum  the  contributions.  As  we  now  show,  these  two 
operations  can  be  carried  out  in  closed  form.  An  important  bonus  of 
having  analytical  expressions  is  the  possibility  to  identify  the 
dominating  matrix  elements  of  A. 

Combining  (C.5)  and  (C.6)  with  (4.6)  implies  that  all  integrals 


(C. 4)  are  combinations  of 


2  Aiv 

/*  dx  xv_1  exp[~  (x  -  J-)]  «  2YV  e  2  K_v(Yy) 


(case  I) 


/“  dx  xV  1  exp[^  (x  +  ^-)]  »  iir  YV  e  2  H^(Yy) 
0  2  x  -v 


for  Im(y)  >  0,  Im(Y  y)  >  0,  any  v 


(case  II) 

(C.8) 


where  K  and  H  denote  the  modified  Bessel  and  Hankel  functions  of  the 
v  v 


first  kind  of  order  v,  respectively.  As  will  become  clear  shortly,  we 


shall  need  only  the  asymptotic  or  nearthe  origin  behavior  the  Bessel 

.  ......  .  .  ..  14 


function,  given  by 


*.  f  \  r  TT-I  1  /2  “  Z 

K  (z)  *  [t-J  e 
-v  2z 

r  »v  JT-. 

„<”(*)  -  e  2  '  4 

-  V  vz 


for  z  -»  » 


(C. 9. a) 


.  '■  '  1 
.  .A  . 


A-8 


k  (2)  -  ~  r( v)  (§) 

-v  2  2 

u(1),  V  1  ,  iuv,z.-v 

H  (z)  - - r(v)  e  (-r) 

~V  IT  2 


for  z  -*  0 


(C.9.b) 


Before  proceeding  further,  we  mention  that  the  conditions  attached  to 
(C.8),  i.e.,  Im(y)  >  0,  etc.  are  always  fulfilled.  Also  because  the  k^- 
integration  in  (C.3)  is  over  the  whole  k^-axis  while  in  (C.8)  it  runs 
over  only  half  the  axis,  we  subdivide  the  k^-integration  into  -®  i 
0  and  0  £  k^  S  ®.  The  respective  contributions  are  denoted  by  the 
and  superscripts  in  (C.10). 

By  virtue  of  (C.8)  and  (C.9),  we  ootain  now  an  estimate  of  the 

magnitude  of  g  (t)  (and  therefore  of  the  t-summed  expression)  by 
*  »  J 

considering  the  "Yy"  argument  in  (C.8).  This  factor  determines  the 

dominating  exponential  or  power  in  the  expression  for  g  (t).  Constant 

*  •  J 

phases,  etc.,  are  omitted  in  (C.10),  but  are  included  in  the  final 
expressions.  Straightforward  algebra  based  on  (C.3)  and  (C.6)  gives 


„±  r  ik2(1)  ,1/2 

■  ['  2TTTT] 


(Case  II) 


Y*  y*  -  2[(2ig  ±  t)(-ik2(1)g)]1/2 


r  2[2g2k2(1)]1/2(1  ±  ||) 


r  2[itgk2(1)]1/2(1  ±  ±Z) 

Y'y*  -  ?[(2ig±t)(ik2(0)g) ]1/2 


for  g  >>  t 


for  g  <<  t 


*  "’■j.  :  ■ .  Vj. y  *  j." 


A-9 


r  2[-2g2k2(0)]1/2(1  +  for  g  >>  t 

-  2[itgk2(0)]1/2  (1  ±  |&)  for  g  «  t 


i  -  i-  ifrrrTTTe^2  <Case  n>  f°r  k  < 

U 

2  2  1/2 

Y*  -  [+  (Case  I)  for  k  >  k  jJl 


4ig  ±  t  ±  v/k. 


Y*y*  -  2[{Dig  ±  t  ±  TT/kQ)(-ig(k^(0)  +  k2£l))]1/2  for  k  <  k^//2 

-  2[i4g2(k2(0)  +  k2(0)]1/2[l  £  ^(t+lt/kG]  for  8  >>  t 

-  2  [g(  t+n/kjJ  i(k2(0)  ♦  k2(l))]U2[l  ±  ~ ]  for  g  «  t 

G 

Y*u*  -  2[(4 ig  ±  t  ±  ir/kG)(Mg(k2(0)  -  k2(1)))]1/2  for  kQ  >  k  /Jl 


2[-4g2(k2(0)  +  k2(1))]1/2[l  +  ^  ( t+ 7r/kQ)  ]  for  g  >>  t 
2[g(  t+Ti/kG)  i( k2(0)  +  k2(  1 ) )  j1  /2[l  +  •  ]  for  g  <<  t 


Y*  -  [ - —k  *£!-*-*■  (--} - ]  (Case  II)  for  k  <  k  //2 

■"V  iTr  *  1  ♦  '/kG]  p 

-  [»  — —  !01  *k  —  — — ]  (Case  I)  fop  k  >  k  /2 

Vofe  1  1  *  */kc1  P 

u 

■  2  (‘sir 1  tJ’/ko,<-  «eib<k2(0)  *  k2(,,)>i’/2  for  ko 2 

G  G 

-  2[(-rJ-  ±  t  +  £-)(•••  — — (k2(0 )  ♦  k2(1)))]1/2  for  kn  >  k//z 
SkG  kG  ^ekr  0  p 

u 


for  all  t 


(C. 10) 


where  k  -  <Wc  and  where  the  square  roots  in  (C.10)  and  throughout  [Eq. 

(2.3)]  are  chosen  so  that  the  real  or  imaginary  parts  are  positive.  We 

have  also  used  in  (C.10),  and  in  the  rest  of  this  Appendix,  eQ-1 ,  and 

the  notation  "Cases  I  and  II"  are  defined  in  (C.8). 

Expressions  (C.10)  yield  a  quick  order-of-magni tude  estimate  for 

g  (t).  Consider  first  the  g  >>  t  regime.  Since  gk(0)  -  gk(l)  -  gk  - 
i » J  P 

6  >>  1 ,  the  elements  belonging  to  classes  a,  6  and  Y  call  for  the 

asymptotic  limit  (C.9.a),  while  those  of  class  6  call  for  the  limit 

(C.9.b).  Therefore,  the  leading  t-independent  terms  in  (C.10)  yield  the 

following  dominating  exponents  in  g  (t)  (Table  1): 

i » J 


class  a  -  e 
class  $  -  e 


2i[2g2k2(1)]1/2 


-2[~2g2k2(0)]1/2 


class  Y 


e2 i[^g2(k2(0)+k2(1))]1/2 

-2[-Hg2(k2(0)  +  k2 ( 1 ) ) ] 
e 


for  k  i  k  /fk 
P 

1/2 

for  k  £  k  //2 
P 


class  6  -  1 


for  g  >>  t  .  (C. 11 ) 

The  symbol  "1"  is  the  class-6  entry  indicates  that  no  exponential 

factors  appear,  only  powers.  Estimates  (C.11)  imply  that  as  B  -*•  ®,  the 

2 

class-a  elements  are  always  the  smallest  [k  (1)  <0,  always].  The 

class-Y  elements  are  larger  than  the  class-a  elements  by  a  factor 

exp[ 2i[ Mg2k2 ( 1 ) ]1 /2k2(0 ) /k2( 1 ) ]  >>  1  for  k  <  k^//!,  and  by  yet  a  larger 

factor  when  k  >  k  //2.  It  is  therefore  justified  to  discard  (set  to 
P 

zero)  the  class-a  elements  for  the  g  >>  t  regime. 

Consider  now  the  complementary  g  <<  t  regime.  Considerations 


similar  to  those  above  yield 


class  a  -  e2lCltgkJ(1)],/2 
class  8  -  e-2Utgk2(0)]’/2 

class  6  -  eMltgO<20>  *  k2(0))]'/2 
_  e-2[itg(k2(0)  ♦  k2(1))]1/2 


for  k  £  k  //2 
P 


for  k  i  k  //2 
P 


class  6  -  1 

for  g  «  t  .  (C. 12) 

Since  the  real  and  Imaginary  parts  of  the  square  root  are  always  chosen 

to  be  positive,  (C.12)  clearly  indicates  that  the  elements  belonging  to 

classes  a,  8  and  Y  are  exponentially  small,  while  the  class-6  elements 

are  of  the  order  1  (i.e.,  fall  off  as  an  inverse  power).  Therefore,  in 

conjunction  with  the  Poisson  sum  (C.3),  there  is  a  sharp  cutoff  in  the 

u-summation  for  the  g1  j(t)  pertaining  to  classes  a,  8  and  Y,  whereas 

for  the  class-6  element  it  is  necessary  to  3um  all  the  way  to  infinity. 

The  cutoff  L  (in  the  u-summation)  occurs  [from  ( C. 10)]  when  2uL/k  -  t  - 

G 

2g  or  L  =  8/ n  >>>  1.  It  is  sharp  in  the  sense  that  over  a  range  AnL  << 

1  the  n-dependent  contributions  diminish  as  e-^". 

The  asymptotic  (8  >>  1)  vanishing  by  the  class-a  elements  implies 

that  in  the  evalution  of  A  1  only  the  eight  elements  of  clases  Y  and  6 
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contribute  .  We  therefore  turn  now  to  the  calculation  of  these  entries 


for  a  fixed  t  and  then  consider  the  summation  mandated  by  the  Poisson 
formula  (C.3).  Straightforward  algebra  gives  for  the  class-Y  elements 


A-12 


g1  2(t)-  P^(g,kli)[|(k2(0)  ♦  k2(  1 )  ]  cos[v^(t)] 

83  j/t)  “  P^8.k(|)[l(k2(0)  *  cos[v^(  t)  ] 

X 

g1  ^(t)  -  P^(g,k|i)[l(3k2(0)  ♦  k2(  1 ) )  ]  i  sln[  v>(  t)  ] 


g3,2(t)  "  p^(g-k|,)  [^k2(°)  +  3k2(  1 ) )  ]  i  sin [ v>( t )  ] 

where  the  labels  >  and  <  refer  to  the  domains  where  k  £  k  //2  or 

P 

kp//2,  respectively,  and 


>,  ,  2C  rk2(0)*k2(  1  i 

P,U.«  )  -  - ; - . 

G 


ti  exp(  in) 


1  /2 


2  ?  ,  A  ,  2 ,  ,  .  1 1  /2 

■g  x  (0)*k  U  )  ] 


1  r  pi  ’2 

«  exp[iTi(-  -  k(l/kG)]  fcxp.-i.-~.-  *  :,♦*  (1))] 


Y>(t)  *  2[-4g2(k2(0)*k2(1  )  )  j 


—t  ♦  n/k  ) 
“  g  o 


?  2  1/2  1/2 
P<  '  i  2iCr  k^(0)  +  k  (U  r  n  exp(  -1  n/2)  1 

^  >  (  8  • K  >  “  K  l "  ij  J  L  2  2  2  1/2  -* 

C  [i<32(k2(0)fck2(1))] 


*  exp[iTr(—  -  k  /kr  , 


. *k2( i ) )  ] ' /2 


v<( t)  »  2i[  kg" 


Note  that  't  (t)  dv  t.a  d-p-r.d  ,.r.  g!  Similarly,  the  class-5  el 
< 


(C. 1 3) 

k  < 


(C. Ik) 

ements 


A- 


g,  .(t)  -  P.(k  )  [2e3TTi/2(l"(g,t)  *  I*(g.t)) 

£ ,  i  o  it  i  i 

*  ^(k2(0)  ♦  k2(l))  ell/2(l“(g,t)  ♦  I*(g.t))] 
g2  3(t)  -  P6(ku)  [2e3^i/2(-i)(I3(g,t)  -  I*(g,t)) 

♦  |  k2(0)  ei,/2(l"(g,t)  -  l|(g,t))] 

g4  ^t)  -  P6 Ck  )  [-21  e3lTi/2(l“(g,t)  -  I^/g.t)) 

♦  |  k2(1)  eilt/2(I~(g,  t)  -  l|(g,t))] 

g4  (t)  -  P6(k  )  2e3iri/2(l'(g,t)  ♦  IjCg.t))  ,  (C.15) 

where 

Pfi(k  )  -  C  ex p [  i (  -  ^  +  nkj|/kG)  ] 

I*(g,t)  ■  [gJT  *  t  +  ir/kG]n  .  (C.  16) 

G 

Expressions  (C.15)  and  (C.16)  are  valid  for  the  whole  k  1  k^  domains. 

The  last  step  is  to  perform  the  t-summation  implied  by  ( C . 3 ) -  As 

pointed  out  before,  for  the  class-6  elements  (C.16),  there  is  a  sharp 

cutoff  at  t  -  2un/kr,  -  2g,  whereas  for  the  class-6  elements  (C.18)  all 

u 

values  of  t  -  2nn/kn  must  be  summed.  The  summations  pertaining  to 

U 

1 8 

( C.  1  3 )  and  (C.15)  are  simply  caried  out  by  using  the  identities 


cos[a(n+1  /2)  ] 


I  (a,Y)  -  l  e  iYW  sin[a(n+1/2)] 
3  n-L 


i 

^  e 


+  nr/2 


r s In [ L ( g+  Y) ]  _  sin[L(g-Y)]  i 
sin(^)  sin(— 


n«-°° 


•  ^  vi  v  -i  Yn  1 

U,  Y)  -  l  e  — — 
1  a+n 


ue 


+  iaY 


3in(ait ) 


I3(a,T)  -  1  —  I1  (a.Tf)  . 
'  da 


Using  ( C.  1  7 )  in  (C.3)  gives  for  the  class-Y  matrix  elements 
> 


A1,2  “  Pa(g’V  [|(k2(0)  +  k2(1))]  Ic(a>,Y) 


A3,4  "  Py(S- k|,>  [^(k2(0)*k2(1))]  Ic(„  t) 
A1,4  “  t^3k2(0  )  +  k2(  1 )  ]i  Is(0>iT) 


A3,2  "  P<r(S*ki|)  [^CK2(0)>3k2(  1 )  Ji  Is(a>,Y)  , 


(C.  17) 


where 


a.  -  i[k  (0)  +  k  (1)]  ir/k.  for  k  <  k  /Jl 

N  UP 

as  -  [-(k2(0)+k2(1  )]1/2ir/kp  for  k  >  k  /Jz 
>  G  p 

and 

Y  “  +2irk(j/kG  . 

For  the  class-6  matrix  elements  we  obtain 
A2  1  “  0  ^of  order  (1/S)]) 

*2,3  ■  5  P6(k  )  e"1"2!  2(.2*y2)(^)  -  k2(0)) 

S.,  -  Wk  >  .-n/2U(.2^)(^)2  -k2(1)| 

^  “  0  [of  order  Qf C  1/6)3 

In  evaluating  (C.20)  we  used  a  «  -1/2  ■*■0(1/3)  »  -0 
I^Ca, Y)  of  Eq.  (C.17),  and  Y  is  defined  in  (C.19). 

vanishes  along  a  curve  wl thin  the 


(C.  18) 

(C.19) 


(C.20) 

5  in  1^ (a,Y)  and 
Note  that  ^  never 
light  cone  k(0)  « 


vanishes  and  A, 


s 


\  • 

a 

w\ 

a 

*V  • 

« 


r 


t 


K" 

%; 


h,  * 

k: 

P 


"■I-'A11  W  t"  LL*  £ 1 11  ’J."  -l1  -1.1 1  l-.1" 1  11  -  ' I'.w?  m;  nr*1?’-11 1  T  "-" 1 '  - 1  •» '  -V  •- ' ■'  *?  T ' "  r  •’ 


A-1  6 


Appendix  D:  Local  Field  Enhancement 

Our  starting  point  is  the  projection  of  the  slow  component  in  the 
field,  E  (X, Z)  [EQ.  (5.1)3,  along  a  vector  v,  which  is  either  x  or  z. 
From  (5.1),  (3-16)  and  (C.20)  one  observes 


v *E  ( x ,z)  -  i  {  l  X.uU)  -  ~ —  [I  X.a .(£)][[  X.a  (£)][£  b  (£)p(£)] 

*  D  £—«  1  N,1  £  1  1  £  £ 

-  7—^—  Cl  X  a,U)][I  b  ( £)p( £)]  } 

2,3  £  5  £ 


X£  ‘  v’p0,*(£)  e 


i[ k^x  +  WQ( £ )z] 


( D .  1 ) 


For  B  >>  1  and  normal  incidence  (k  -  0),  the  entries  to  (D.1)  are  [Eqs. 
(2.6),  (3-3),  (k  1  6)  and  (C.2)] 


lid) 


at(£) 


a^(£) 


b2(£) 


-^7=  -J-  expCC-|l|/2B3 

e  /[IT 

lkG  exp(-iir/ii)  |£|  exp[-|£|/2g] 
kG  expL-in/^]  £  ex  pC  - 1  £  1 2 

C  ikQ  |  £|  (-1  )l  exp[-[£j/2B] 


£  |  >>  1 


b4(£) 


C  kQ£ ( -1 ) £  ex p[ - J £ | /2B] 


( D.2 ) 


1  M 

and  from  (2.3)  one  deduces 

u(0)  =*  ♦  i  exp[-i/gk(0)  +  n/ii]  - k- — — — ~v/p  — ~  •  (D.3) 

[k(0)  +  k(1)]J  /2^g 

We  now  specialize  (D.1)  to  points  T  and  B  in  the  notation  of  Fig. 
1.  Using  the  definitions  (2.3),  one  obtains 
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or  explicitly 


2 


for  r  <<  1 


M  (t)  | 2  -LU2i  | - !LXJ _ | 

3  1  D2  1  (/r\  /T)3/2' 


■0  1 


1_  k_(0_)  2 

rv2  2  it  S  T.x 


for  r  >>  1 


2  2 

i:if 


Iz-E  (T)|2  .  ’ - .".K  (°> _  XT,? 

1  s'  Jl  „2  _2r  1  .2,2  2ttS 


rv4-  r  1  t-ii 

D  r  -  r  ] 


1 2  1  4k  (0 ) 


|2‘Es(B)l  d2  2tt B  p2(i:r2)2  *2  *B,* 


2  2 
n  r 


where 


“  l  /PT  M>*  exp[-|£|/6]  -  -1 

£-- *  /g 


T ,  x 


[r] 

I 

|£|>1  /jTf 

[r] 


—  exp[-| £|/2Bl  ~  /r  , 


IT,z  *  .  8,2  exPt"l  1|/26]  -  r  3 


B ,  x 


£»-[r] 

I  /nT(-l)£  exp[-|l|/2B]  -- 


£,*  -  00 


/b 


1  <<  r  <  2B 


B ,  z 


I  r  ( ~1 ) 4  exp[-|£j/2B]  -  B 


£«-oo 


( D.7 ) 


(D.8) 


We  added  in  (D.7)  an  order-of-magni tude  estimate  for  the  various  sums. 
Using  (D.6),  we  obtain  (5.3) • 
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Schematic  display  of  the  sinusoidal  [z  -  g  cos(2nx/d)]  grating  and 


the  notation  used  in  this  work.  The  incident  light  (En  and  the 

Op 

energy  Bragg  reflection  are  indicated.  The  points  nT"  and  "B"  are 
those  for  which  the  local  field  enhancement  is  evaluated. 

Display  of  the  surface-plasmon  flat  bands,  Eq.  (^.^).  The  light 


cone  is  indicated  left  of  the  diagonal  line,  and  k^ 
volume  plasma  wave  vector. 


u  ,  is  the 
p/c 


\\  *. 
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